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Abstract 

S 

Starting from the hypothesis of scaling solutions, the general exact form of 
the scalar field potential is found. In the case of two fluids, it turns out to be 
a negative power of hyperbolic sine. In the case of three fluids the analytic 
^ ■ form is not found, but is obtained by quadratures, 

m ! PACS: 98.80.Cq, 98.80.Hw, 04.20.Jb 
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Introduction 

In cosmological theories containing scalar fields with a self- interaction potential V(<p), it is 
sometimes possible to reconstruct the required scalar field potential for a simple cosmological 
(3jT)! solution. In the context of inflationary theory, this approach was used by various authors 
0, §, §, [|], [|, § [0] § §, 0, and was primarily interested in the behaviour of 
solutions containing only scalar fields undergoing inflation. The more recent invocation of a 
scalar field as a dark matter source responsible for accelerating the universe today, under the 
pseudonym of "quintessence", fLl] , |13|| , flU] , is mathematically almost identical, 



X 



but places different requirements on the solutions. In particular, it is of interest to find 
solutions which contain both perfect fluids and scalar fields. In this new context, interesting 
applications of the reconstruction approach were made by Chiba and Nakamura [TE\ and 
Saini et al. fF< 



In this paper, we seek exact cosmological solutions for a universe containing a 
perfect fluid and a scalar field. We start from the assumption that the energy density of the 
scalar field scales as an exact power of the scale factor: p^ = K a~ n , which is equivalent to 
imposing an equation of state linking the pressure and density, of the form p v = wp^ with 
constant w — n/3 — 1. For a flat Friedmann universe, it is then possible to find an explicit 
exact form of the potential in terms of n, H Q , and fl m o- The form of the resulting general 
solution has instructive features which will be discussed below. 
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Derivation of the Potential 



Consider a cosmological fluid with two non-interacting components: perfect-fluid matter 
and a scalar field ip with potential V((p). In the flat universe case, we have the equations 

3H 2 = g( Pm + Pip ) (i) 

(p + 3H<p + V'(ip) = (2) 

= \p 2 + V(<p) (3) 
Pm = Da~ m , (4) 

where H = a/a, where a(t) is the expansion scale factor is the Hubble expansion rate, and 
overdot denotes differentiation with respect to the comoving proper time t; Q — 8ttG/c 2 and 
V = dV/dip. The constant m depends on the type of perfect fluid present. 

We normalize the present value of the scale factor to a — 1, without loss of generality, 
and for brevity denote the present matter-density parameter, f2 m0 , by Q , and define 

D = 3Hfo /G . (5) 

If we make the assumption that 

p^ = Ka~ n , (6) 
with n < m, so the scalar field can dominate at late times, and define 

k = 3H 2 (i - n )/g , (7) 

then from Eqs. (2), (3) and (6) we obtain 

f = T«" n . (8) 



Since 



we have 



^ = Ha^ , (9) 
dt da' w 



dip \ 2 K s/n 



da J ' 3H 2 Q a n - m + 2 + (1 - Q )a 2 



(10) 



which gives 
ip(a) = J 



1 — Vtn \fnda 



Q s/VL a n - m+2 + (1 - Q )a 
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2y" . i / / 1 — m ~ n i 

-arcsmh W — a 2 + tp in . (11) 



\/G(m — n) IV ^0 
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Returning to the potential, we get from Eqs. (3) and (9) 

V = Ko* -\? = ^{l- Qo) (l - I) a- • (12) 
The scale factor can be easily eliminated, giving eventually 



^) = Y (1 " fio) ( X 



v 'V 6M O 
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smhiVG—^=-((p-<f in ))) . (13) 
This expression, in the case of dust (m = 3), coincides with the one presented by (18 



(which is in turn a particular case of the treatment of |]19|)as well as with the one by |20] 



A general discussion about exact solutions for Friedmann equations, which includes ours as 



a particular case, can be found also in J7J, ||, ||. In |T8| and [1{J, one can find explicit 



solutions for a(t) and (f(t), as well as an extensive discussion. We remark only that this form 



of potential is good for a tracker solution [13], [14|. Indeed, straightforward computation of 



the function T = V"V/(V') 2 , introduced and discussed in these papers, gives 

m — n f f m — n \\ 2 , . 

T = 1 + — sech — > 1 (14) 



2n \ \ \/n 
as required by the tracking condition 



Our derivation differs from flip] , |T9|1 and [£(J because it is simpler and is generalised to 
include all perfect fluid equations of state (other than the n = 6 case which would correspond 
to a pure scalar field with no potential). Moreover, it proves that this form of V(<p) is the 
unique solution, if condition (6) is imposed. 



Discussion 

The first interesting feature of Eq. (13) is that the slope (as well as the amplitude) of 
V(ip) depends on n. This means that it is impossible to obtain a scaling solution, with the 
same potential slope, when passing from a radiation-dominated (m = 4) epoch to a matter- 
dominated (m = 3) epoch. Even if we assume that n changes in such a way that the slope 
remains constant, the coefficient (1 — n/Q) changes. Moreover, the effective equation of state 
of the scalar field also changes and there is no physical mechanism for it to be influenced in 
this way by the dominating type of matter. Although this situation seems to be unphysical, 
it is merely an artifact of having sought a solution containing a single perfect fluid. The full 
solution must be found by including dust, radiation and scalar field from the outset, not by 
joining the radiation + scalar solution to the dust + scalar solution. If this is done for the 
case of dust + radiation + scalar field, with the same arguments as before, it is easy to 
derive 



J V Q ^tt r a n - 2 + IV™- 1 + %a 2 Vm 1 ' 

V = ?Sn v (l - a"» , (16) 



g * v e 
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where Q r , Qd and Vt^ are the present values of the radiation, dust and scalar field density 
parameters, respectively, so Q, r + fid + Q v = 1. These equations give a parametric represen- 
tation of V(ip), which cannot be solved analytically with simple functions, but is otherwise 
perfectly well defined and interpolates between the dust and radiation solutions of Eq. (13). 
In our universe today, Q r << Qd, but at early times the radiation term Q r a n ~ 2 dominates 
the dust and cannot be dropped. 

A numerical example illustrates the situation: let us set n = 1, Q v = 0.7, Q r = 0.0001, 
Qd — 0.3 — Q r , 3Hq — 1, Q — 1. Fig. 1 shows that on the first part the "true" potential fits 
well with Eq. (13) and m = 4, while Fig. 2 shows that in the late regime the fit should be 
done with m = 3. 

The simple and attractive form of Eq. (13) is lost, although it might be recovered by a 
suitable choice of the exponent: via some weighted mean with of Q r and Qd- But there is 
no way of doing this other than a fit of the numerical values. Moreover, there is no reason 
why the exponent should depend on the particular values of Q r and Qd- Also, since in this 
case the scaling feature of the solution is only approximate, the tracker behavior could be 
affected. 

Simple numerical evaluations (with the parameters within the allowed range) show 
that the "true" potential is very well approximated by 

V=£-7, (17) 

with a, j3, 7 positive constants depending on the parameters; in particular, (3 2n/{m — n). 
The negative additive term, 7, is an artifact of the approximation over a finite range of (p. 
This result, already found by |l^] for the dust case, shows that the inverse-power potential 
is effectively equivalent to Eq. (13), so that it is possible to apply to this situation all 
the known results about tracker solutions. On this point, it is interesting to note that 
this approximation is very good over the whole range of a from zero to the present-day 
value (ao = 1 according to our normalization). The asymptotic exponential behavior of the 
hyperbolic potential is therefore important only in the very far future and does not affect 
the dominance of scalar field now or its behaviour in the recent past. 

As a final remark, we should stress that our arguments are based on the arbitrary 
assumption of Eq. (6). We have shown elsewhere |^T| that other forms of exponential 



potential are perfectly able to reproduce observational data, but of course in these cases w 
is not constant. It is interesting also to note that, in one of the cases treated in that paper, 
w is almost perfectly constant, and yet the form of the potential is substantially different 
from that of Eq. (13). 
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Figure Captions 

Figure 1: Early-time regime: the dashed curve is the n = 1 scalar field potential with 
dust (m = 3); the solid curve is the n = 1 scalar field potential with radiation (m = 4) and 
lies closer to the points plotting the full numerical solution for dust plus radiation and n = 
1 scalar field potential. 

Figure 2: Late-time regime: the dashed curve is the n = 1 scalar field potential with 
radiation (m = 4); the solid curve is the n = 1 scalar field potential with dust (m = 3) and 
lies closer to the points plotting the full solution for dust plus radiation and n = 1 scalar 
field potential. 
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Fig. 1 
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